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The low energy behavior of a recently proposed model for the massive analytic running coupling of 
QCD is studied. This running coupling has no unphysical singularities, and in the absence of masses 
displays infrared enhancement. The inclusion of the effects due to the mass of the lightest hadron 
is accomplished by employing the dispersion relation for the Adler D function. The presence of the 
nonvanishing pion mass tames the aforementioned enhancement, giving rise to a finite value for the 
running coupling at the origin. In addition, the effective charge acquires a "plateau-like" behavior 
in the low energy region of the timelike domain. This plateau is found to be in agreement with a 
number of phenomenological models for the strong running coupling. The developed invariant charge 
is applied in the processing of experimental data on the inclusive r lepton decay. The effects due to 
the pion mass play an essential role here as well, affecting the value of the QCD scale parameter A 
extracted from these data. Finally, the massive analytic running coupling is compared with the 
effective coupling arising from the study of Schwinger-Dyson equations, whose infrared finiteness 
is due to a dynamically generated gluon mass. A qualitative picture of the possible impact of the 
former coupling on the chiral symmetry breaking is presented. 
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I. INTRODUCTION 

The theoretical analysis of strong interaction processes 
at low energies represents a long-standing challenge for 
Quantum Chromodynamics (QCD). Whereas the discov- 
ery of asymptotic freedom [lj was followed by the rapid 
development of perturbative tools for the detailed study 
of the ultraviolet region, a reliable method for descrip- 
tion of hadron dynamics in the infrared domain is still 
missing. Given that many important QCD phenomena, 
such as hadronization, quark confinement, chiral sym- 
metry breaking, and dynamical mass generation, are in- 
frared in origin, one resorts to the variety of models, in an 
attempt to obtain a consistent quantitative description of 
the low energy dynamics. 

The renormalization group (RG) method @, Q plays 
a fundamental role in the framework of Quantum Field 
Theory (QFT) and its applications. In the case of QCD, 
in order to describe the physics in the asymptotical ul- 
traviolet region, one basically applies the RG method to- 
gether with perturbative calculations. In this case, owing 
to the asymptotic freedom, a priori unknown RG func- 
tions can be parameterized by power series in the strong 
running coupling. Eventually, this leads to approximate 
solutions of the RG equations, which are used in the 
quantitative analysis of the high-energy processes. How- 
ever, such solutions possess unphysical singularities in 
the infrared domain, contradicting the general principles 
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of the local QFT, and complicating the theoretical de- 
scription and interpretation of the intermediate- and low- 
energy experimental data. Nonetheless, these difficulties, 
being artifacts of the perturbative treatment of the RG 
method, can be circumvented by judiciously incorporat- 
ing nonperturbative information about the hadron dy- 
namics at low energies. 

It is worth mentioning several well-known examples 
of such "synthesis" . The short-range part of the static 
quark-antiquark potential can be calculated perturba- 
tively , while its linear confining behavior at large dis- 
tances is corroborated by both the lattice results (see 
recent papers 0) and the string hadron models (see, 
e.g., book and references therein). These two in- 
puts complement each other and form the so-called "V- 
scheme" Q for the QCD effective charge, which has 
proved to be successful in describing hadrons as bound 
states of quarks 8] . The so-called "I-scheme" [t| is con- 
structed along the same lines. Here, the perturbative 
results are supplied with the large distance behavior of 
the running coupling, coming from the lattice study 0] 
of the topological structure of the QCD vacuum. In- 
terestingly enough, both aforementioned schemes, al- 
though being based on different assumptions, predict 
a similar infrared behavior for the strong running cou- 
pling. Furthermore, the latter also agrees with that of 
the model for the QCD analytic invariant charge devel- 
oped in Refs. 0, ^3 ( see a l so Refs. fl3L Il4| for the de- 
tails). There is also a number of methods which proceed 
from the general properties of the perturbative power 
series for the QCD observables in the framework of the 
renormalization group formalism. For example, these are 
the "optimized perturbation theory" 0,0], the method 



2 



of effective charges Il7|. the Brodsky-Lepage-Mackenzie 
convergence criterion [18j, the "optimal conformal map- 
ping" method ^ , and the RG improvement of pertur- 
bative calculations [20L l2l| . 

Another important source of nonperturbative informa- 
tion is provided by the relevant dispersion relations. The 
latter, being based on the "first principles" of the the- 
ory, supply one with the definite analytic properties with 
respect to a given kinematic variable of a physical quan- 
tity in hand. The idea of employing this information 
together with perturbative treatment of the renormal- 
ization group method forms the underlying concept of 
the so-called "analytic approach" to QFT. It was first 
proposed in the framework of Quantum Electrodynam- 
ics (QED) and applied to the study of the invariant 
charge of the theory |23 • Here the principle of causality 
implies the Kallen-Lehmann spectral representation for 
the QED running coupling. Hence, the latter has to be an 
analytic function in the complex g 2 -plane with the only 
cut along the negative 1 semiaxis of real q 2 . A number 
of authors (see, e.g., Ref. [23J) have argued that a simi- 
lar method can also be useful for studying non-Abelian 
theories. Eventually, procee ding from these motivations, 
the "dispersive approach" |24[ and the "analytic ap- 
proach" |25j to QCD have been developed. According to 
the former one, the nonperturbative effects of the strong 
interaction can be reliably captured at an inclusive level 
by means of a quantity, which constitutes the effective 
extension of the perturbative running coupling to the 
low energy scales. The analytic approach to QCD has 
been successfully applied to the study of the strong run- 
ning coupli ng [ la. l2q . perturbative series for the QCD 
observables |26(, and some intrinsically nonperturbative 
aspects of the strong interaction [ll|, [lj, |23 . Some of the 
main advantages of the latter approach are the absence 
of unphysical singularities and a fairly good higher-loop 
and scheme stability of the outcoming results. Besides, in 
the framework of the analytic approach the continuation 
of the "spacelike" theoretical predictions for the QCD 
observables into the timelike domain, that is crucial for 
handling the relevant experimental data, can be carried 
out in a self-consistent way pgj . 

In general, the effects due to the masses of light 
hadrons (such as tt meson) can be safely neglected only 
when one studies the strong interaction processes at large 
momenta transferred. For example, in order to relate the 
perturbative results with the high energy experimental 
data on the electron-positron annihilation into hadrons, 
the massless approximation of the dispersion relation for 
the Adler D function may be used (see Section [n] for the 
details). But for the hadron dynamics in the infrared do- 
main the mass effects become substantial. Apparently, 
this is important for the description of the low energy 



experimental data on the inclusive r lepton decay. Both, 
the results of perturbation theory and the dispersion rela- 
tion for the Adler D function with the nonvanishing pion 
mass, are vital here for properly processing these data. 
However, no such mass effects have been taken into ac- 
count within the analytic approach to QCD so far. 

The primary objective of this paper is to include the 
effects due to the pion mass into the analytic approach 
to QCD. The incorporation of such mass effects is stud- 
ied on the example of the model for the analytic run- 
ning coupling developed in Refs. 0,0]. Therein, the 
imposition of the analyticity requirement has eventually 
resulted in the infrared enhancement (i.e., the singular 
behavior at q 2 — 0) of the invariant charge in hand. In 
general, one might anticipate that the presence of masses 
affects the low energy behavior of the strong running cou- 
pling. Indeed, as we shall see, the aforementioned singu- 
larity is tamed down by the pion mass, thus giving rise 
to a finite infrared limiting value for the QCD effective 
charge. Apparently, it is important to apply the devel- 
oped model to the description of those sets of experimen- 
tal data, which display a particular sensitivity to the in- 
frared behavior of the QCD running coupling. It is also 
of significant interest to examine, even at a qualitative 
level, the applicability of the obtained invariant charge 
to the study of the chiral symmetry breaking through 
Schwinger-Dyson equations. 

The layout of the paper is as follows. Section [H] is 
devoted to the description of the strong interaction pro- 
cesses in spacelike and timelike domains. This material 
sets up the stage for the subsequent analysis of the mass- 
less and massive cases. In Section IIIII the analytic ap- 
proach to QCD is overviewed, with a particular emphasis 
on the massless model for the invariant charge of |Iull2| . 
The effects due to the pion mass are incorporated into the 
latter approach in Section llVl The basic features of the 
massive strong running coupling in spacelike and time- 
like regions are also studied therein. In Section [V] the 
developed model for the invariant charge is applied to 
processing the experimental data on the inclusive r lep- 
ton decay, a reasonable estimation of the QCD scale pa- 
rameter A being obtained. In Section IVII the derived 
massive analytic charge is compared with the effective 
charge arising from the study of the Schwinger-Dyson 
equations, whose infrared finiteness is due to a dynam- 
ically generated gluon mass [2^. A qualitative picture 
of the possible impact of the former charge on the chiral 
symmetry breaking is presented. In Conclusions (Sec- 
tion IVII|) the basic results are summarized and further 
studies within this approach are outlined. 



II. STRONG RUNNING COUPLING IN 
SPACELIKE AND TIMELIKE REGIONS 



1 A metric with signature (—1, 1, 1, 1) is used, so that positive <y 2 
corresponds to a spacelike momentum transfer. 



The consistent description of hadron dynamics in time- 
like (Minkowskian) and spacelike (Euclidean) regions re- 
mains the subject of intense studies. The strong inter- 
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action processes involving the large spacelike momentum 
transfer q 2 > (for instance, the deep inelastic lepton- 
hadron scattering) can be examined perturbatively in the 
framework of the RG method (see, e.g., Ref. [3(j). How- 
ever, in order to handle the processes which depend on 
the timelike kinematic variable s = — q 2 > (for ex- 
ample, hadronic width of the r lepton decay or total 
cross-section of the electron-positron annihilation into 
hadrons), one first has to relate the results of perturba- 
tion theory with the measured quantities. Obviously, the 
question what is the expansion parameter for the QCD 
timelike processes arises at this stage [3l|. 

An indispensable method for the analysis of the 
strong interaction processes in the timelike domain has 
been pro posed by Adler [^2, and further developed in 
Refs. |33|,|3_^]. In particular, it was argued that the log- 
arithmic derivative of the hadronic vacuum polarization 
function II (q 2 ) 



D(q 2 



dll(q 2 ) 
dlnq 2 



(1) 



which is also known as the Adler D function, provides 
a firm ground for comparing the perturbative results 
with the experimental data on the e + e~ annihilation into 
hadrons. Specifically, the dispersion relation |32( 

R(s) 



D(q 2 ) = q 2 



- ds 



(2) 



embodies the required link between the measurable ratio 
of two cross-sections |33 



R(s) 



a (e + e — > hadrons; s) 
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and the Adler D function, which can be calculated per- 
turbatively. In Eq. s denotes the center-of-mass en- 
ergy of the annihilation process. Thus, one can continue 
the perturbative results for D(q 2 ) into the timelike do- 
main by making use of the relation inverse to Eq. @ 
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where the integration path lies in the region of analyticity 
of the function £>(—£), see also Refs. [33l l3fil | . 

So far, there is no systematic method for calculating 
the Adler D function. Nevertheless, its asymptotic ul- 
traviolet behavior at q 2 — > oo can be computed pertur- 
batively. There, the effects due to the masses of light 
hadrons can be neglected, and the Adler D function of 
Eq. Q is usually approximated by the power series in 
the strong running coupling a s (q 2 ) 



D(q 2 ) = N C J2Q} [1 + d(Q 2 )] , 
/ 



(5) 



where N c = 3 is the number of colors, Qf stands for the 
charge of the /-th quark, 



d(q 2 ) ~ di 



[«s(<Z 2 )" 
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\a s (q 2 )] 
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d\ = 1, g?2 — 1-9857 — 0.1153rtf, and rif is the number of 
active quarks, see Refs. |3JJ |37j for the details. 

Thus, in order to compare the perturbative results with 
the timelike experimental data, one first has to perform 
on Eq. © the integral transformation given in Eq. 
It is worthwhile to underscore that this procedure dis- 
torts the perturbative power series for the Adler D func- 
tion drastically, since both real and imaginary parts of 
the running coupling a s (q 2 ) contribute to Eq. Ulti- 
mately, the continuation presented in Eq. |0J results in 
a "non-power" expansion for i?(s), and even in the deep 
ultraviolet asymptotic \q 2 \ — > oo the functions D(q 2 ) and 
R(s) are different, starting from the three-loop level, due 
to the so-called 7r 2 -terms. Nonetheless, the "naive" ex- 
trapolation of the strong running coupling to the timclikc 
domain a(s) = a s (\q 2 \) is also allowed for the perturba- 
tive expansion of Eq. JBJ, but only if one restricts oneself 
to the deep ultraviolet limit \q 2 \ — > oo of the one- or two- 
loop levels (see Refs. [H Hi 1^1 H H |M HJ for the 
details). 

Since the integral transformation Q of the perturba- 
tive results has to be carried out every time one deals 
with the timelike strong interaction processes, for prac- 
tical purposes it is convenient to define [2^| the timelike 
effective charge S(s) in the same way, as R(s) relates 
with D(q 2 ): 



a(s) = lim 

2ni s^o + 



s-\-ie 



( r\ d < 
&(-() —■ 



(7) 



In what follows the strong running coupling in the space- 
like domain is denoted by a(q 2 ), and in the timelike do- 
main by a(s). Obviously, the inverse relation between 
these effective charges 2 



a(q 2 ) 



a(s) 



4ml (s + q 2 )'- 



■ ds 



(8) 



holds as well 3 . It is important to emphasize that for a 
detailed description of the infrared hadron dynamics the 
pion mass cannot be neglected in Eqs. @ and JHJ. 

Apparently, for the self-consistency of the method de- 
scribed above, one first has to bring the perturbative ap- 
proximation for the Adler D function in Eq. (jSJ to con- 
formity with the dispersion relation of Eq. @. This is of 



2 The ca se of the massless pion = was studied in Refs. l2St 

MM 

3 The relations J7J and JSJ are not valid for the perturbative run- 
ning coupling oi s (q 2 ) because of the unphysical singularities of 



the latter, see Section lllll for the details. 
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a great significance when one intends to study the QCD 
experimental data in the intermediate- and low-energy 
regions. Indeed, the integral representation in Eq. (J2J 
implies the definite analytic properties in q 2 variable for 
the Adler D function. For example, in the masslcss limit 
(to^ = 0), it has to be an analytic function in the com- 
plex (/ 2 -plane with the only cut — oo < q 2 < along the 
negative semiaxis of real q 2 . However, the approxima- 
tion of the right hand-side of Eq. JSJ by the perturba- 
tive expansion in the strong running coupling given in 
Eq. © obviously violates this condition. Nevertheless, 
this discrepancy can be eliminated in the framework of 
the analytic approach to QCD, which is discussed in the 
next section. 



III. MASSLESS ANALYTIC RUNNING 
COUPLING 

As has already been mentioned in the Introduction, the 
dispersion relations play a central role in the description 
of hadron dynamics. Indeed, the general principles of the 
local QFT (such as causality, spectrality, unitarity) are 
captured by the relevant integral representations. These 
are, for instance, the dispersion relation for the Adler 
D function and the Jost-Lehmann-Dyson represen- 
tation ^tj f° r the structure function of the deep inelastic 
lepton-hadron scattering processes. In turn, the disper- 
sion relations provide one with a certain nonperturba- 
tive information about the quantity in hand, in particu- 
lar, with the definite analytic properties in the kinematic 
variable. Undoubtedly, the latter should be taken into 
account when one intends to venture beyond the realm 
of perturbation theory. 

It has recently been argued pil l25l | that for the QCD 
invariant charge a(q 2 ) the Kallen-Lehmann spectral rep- 
resentation 



a{q 2 ) = I J^ dn 



(9) 



must hold in the absence of masses. The condition @ 
is identical to that needed 4 for bringing the perturba- 
tive approximation of the Adler D function in Eq. (jSJ) 
to conformity with its dispersion relation (J2J), also en- 
forcing the validity of Eqs. and JSJ|. However, there 
are several ways to incorporate the analyticity require- 
ment of Eq. © for the QCD running coupling into the 
RG formalism. In other words, the perturbative asymp- 
totic behavior of a s (q 2 ) when q 2 — > oo, together with 
the integral representation Q , is not enough to uniquely 
determine the relevant spectral density g(u). Eventu- 
ally, this ambiguity has given rise to different models for 
the strong running coupling within the analytic approach 



to QCD (discussion of this issue can also be found in 

Refs. miii mm si). 

This section is devoted to a brief overview of one 
of the massless models for the QCD analytic invariant 
charge 0, 0] • This model shares all the advantages of 
the analytic approach, namely, it contains no unphysical 
singularities, and displays good higher loop convergence 
and mild dependence on the subtraction scheme. Besides, 
the running coupling of Refs. 0,0] was successful in the 
description of a wide range of QCD phenomena 0H3. 
Furthermore, it is of a particular interest to note that 
this model has recently been re-derived, proceeding from 
completely different motivations |l3| . 

In the framework of perturbation theory the RG equa- 
tion for the QCD invariant charge a(fi 2 ) = g 2 (/j, 2 )/ (4ir) 
at the £-\oop level takes the form 



rima^V) 
d In fj, 2 



e-i 

3=0 



47T 



(10) 



Here ai^ip?) denotes the ^-loop perturbative running 
coupling, [3j stands for the (3 function expansion coeffi- 
cient (/3 = 11 - 2n f /3, /3i = 102 - 38n { /3, ... ), and n f 
is the number of active quarks. It is well-known that 
the solutions to Eq. (|10f> have unphysical singularities in 
the infrared domain at any loop level. Specifically, the 
Landau pole appears at the one-loop level, whereas the 
higher loop corrections introduce additional singularities 
of the cut type into expression for the QCD invariant 
charge. In turn, this contradicts the fundamental princi- 
ples of the local QFT, violating the representation given 
in Eq. ©. 

In order to resolve this difficulty, in the framework of 
the developed model [TTL IT^ the analyticity requirement 
was imposed on the (3 function perturbative expansion 5 



dlnq&V) 
d ln/x 2 




4tt 
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In this equation ain(^ 2 ) is the £-\oop analytic invariant 
charge and the braces { . . . } denote the "analytiza- 
tion" of the expression contained in them |2^| : 



U(q 2 )} 
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-A(-a + is) 



A(-a - ie) 
da 



<J + q* 
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In the limit of the masslcss pion m„ 



5 Unlike the Shirkov— Solovtsov model |25| . where the analyticity 
requirement 1121 was imposed on the perturbative running cou- 
pling a s (q 2 ) itself. In turn, this has led to a spectral density 
somewhat different from that of Eq. 1171 . and, consequently, to 
different properties of t he QCD effect ive charge in the infrared 
domain, see, e.g., Refs. for the details. 



It is worth noting here that the way of incorporating 
the analyticity requirement into the RG method given in 
Eq. 1)11(1 is consistent with the gene ral definition of the 
QCD invariant charge, see Refs. [l4L l41| . 

At the one-loop level the RG equation (|llfl for the 
analytic invariant charge can be solved explicitly |ll| : 



47T Z - 1 



Po z In z 



Si 
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At the higher loop levels only the integral representation 
for the analytic running coupling has been derived. So, 
at the l-loop level the solution to Eq. ijTTJ acquires the 
form mEI: 
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exp 



V W (a) lnfl 



(3q z In z 

where V&{o) = U^{a) - TZ^(a) and 
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3=0 v ' 

- «f)(-a + K )f +1 }. 
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The obtained massless running coupling 1(14(1 has the 
correct analytic properties in the q 2 variable demanded 
in Eq. (JSJ, namely, it has the only cut q 2 < along 
the negative semiaxis of real q 2 . In particular, the latter 
follows from the Kallen-Lehmann integral representation 
that holds for the invariant charge ((141) : 
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In this equation pW(a) denotes the £-loop spectral den- 
sity 
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(18) 



(19) 



is the one-loop spectral density. In the exponent of 
Eq. HI 7(1 the prin cipal value of the integral is assumed 
(see Refs. [Hill for the details). 

The massless analytic running coupling of Eq. 1(14(1 pos- 
sesses a number of appealing features. First of all, it has 




FIG. 1: The integration contour in Eq. £J for the massless 
case. The physical cut of the strong running coupling ct(—Q 
(see Eq. ©) is shown along the positive semiaxis of real £. 



no unphysical singularities at any loop level, and con- 
tains no adjustable parameters 6 . Thus, similarly to the 
perturbative approach, the QCD scale parameter A re- 
mains the basic characterizing quantity of the theory. In 
addition, the invariant charge ((14(1 incorporates the ul- 
traviolet asymptotic freedom with the infrared enhance- 
ment in a single expression, which plays an essential role 
in applications of the developed model to the descrip- 
tion of the quenched lattice simulation data (18l l27L |44j . 
Moreover, this analytic running coupling has universal 
asymptotics both in the ultraviolet and infrared regions 
at any loop level, and displays a good higher loop and 
scheme stability. The detailed analysis of the properties 
of the invariant c harg e (11411 and its applications can be 
found in Refs. Hi HH Eg . 

As has been discussed in Section |n] for the consistent 
description of a number of strong interaction processes 
one has to employ the continuation of the QCD effective 
charge to the timelike region, in the way given in Eq. 10. 
For the massless case under consideration it is convenient 
to choose the integration contour in Eq. Q in the form 
presented in Figure Eventually, this leads to the fol- 
lowing extension of the invariant charge of Eq. 1(141) to 
the timelike domain 1121 
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6 It is wort h noting here that the Shirkov-Solovtsov running cou- 
pling [2.5ij has no adjustable parameters, either. So, both these 
models are the "minimal" ones in this sense. 
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FIG. 2: The one- loop massless analytic running coupling 
in the spacelike (Eq. (tTol . q 2 > 0) and timelike (Eq. 
s = — q 2 > 0) regions. The values of parameters are: 
A = 508 MeV, n f = 2 active quarks. 



sions have identical behavior determined by the asymp- 
totic freedom. However, there is an asymmetry between 
them in the intermediate and low energy regions. Thus, 
the relative difference between these effective charges is 
about several percent at the scale of the Z boson mass, 
and increases when approaching the infrared domain. 
Evidently, this circumstance has to be taken into ac- 
count when one handles the experimental data (see also 
review [l4| and references therein for the details). 

It is worthwhile to emphasize that the mass effects have 
not been included in the formulation of the models for the 
strong running coupling in the framework of the analytic 
approach to QCD, so far. Thus, the obtained results 
can be applied, for example, to the study of the exper- 
imental data at high energies, where the masses of the 
lightest hadrons can be neglected, the pure gluodynam- 
ics, and the quenched lattice simulation data (see also 
Refs. 0,0). However, for the detailed description of 
the infrared hadron dynamics, the mass effects have to 
be incorporated into the analytic approach to QCD. The 
next section is devoted to this task. 



where s = —q 2 > 0, and the spectral density pW(cr) is 
defined in Eq. Ijf 7|l . The obtained result supports the hy- 
pothesis due to Schwinger |4o.l47| concerning the propor- 
tionality between the (3 function and the relevant spectral 
density (see also Ref. p8|l. 

The one-loop effective charge of Eq. (|2U[l has the fol- 
lowing asymptotic in the high energy limit s — *■ oo: 
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(21) 

On the one hand, this running coupling has the correct 
ultraviolet behavior, determined by the asymptotic free- 
dom. On the other hand, the so-called 7r 2 -terms have 
also appeared in the expansion (j21|) . As it was noticed in 
Section [Hj these terms play a key role in the description 
of the strong interaction processes in the timelike domain. 
It is interesting to note that, similarly to the "spacelike" 
running coupling in the massless case of Eq. (|13|l . the 
one-loop effective charge (|20|) also has an enhancement 
in the infrared domain (see Refs. [j^ll4] |): 
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(22) 



However, the type of this singularity differs from that 
of the invariant charge in Eq. (|13|1 by the logarithmic 
factor. Nevertheless, it is precisely this feature of the 
timelike running coupling that enables one to handle 
the integrals of a specific form over the infrared region, 
and in particular, to process the experimental data on 
the inclusive semileptonic branching ratio for the case 
of the massless pion. In turn, the latter provides one 
with the relevant estimation of the QCD scale parameter 
A = (508 ± 6f ) MeV, see SectionEJfor the details. 

The plots of the functions a^n(q 2 ) and ain (s) are 
shown in Figure [21 In the ultraviolet limit these expres- 



IV. MASSIVE ANALYTIC EFFECTIVE 
CHARGE 

As has been noticed in previous sections, the tt meson 
plays a crucial part in the description of the strong inter- 
action processes at low energies. So far the main thrust 
of the analytic approach to QCD has focused on elimi- 
nating intrinsic difficulties of perturbation theory, such 
as the unphysical singularities of the strong running cou- 
pling (see Section llTT|) . On the other hand, mass effects 
within this formalism remain largely unexplored, thus 
far. Therefore, the objective of this section is to incor- 
porate the effects due to the pion mass into the analytic 
approach to QCD. 

Evidently, the original dispersion relation for the Adler 
D function [32] (see Eq. (J2J) with the nonvanishing mass 
of the 7r meson is the proper object to study here. Indeed, 
Eq. J2| implies definite analytic properties in the q 2 vari- 
able for D(q 2 ). Namely, the latter has to be an analytic 
function in the complex g 2 -plane with the only cut begin- 
ning at the two-pion threshold — oo < q 2 < —Am 2 along 
the negative semiaxis of real q 2 . However, its approxima- 
tion in Eq. © violates this condition due to the spurious 
singularities of the perturbative running coupling a s (q 2 ). 
Nevertheless, this disagreement can be avoided by im- 
posing the analyticity requirement of the form 7 



x(a) 



da 



(23) 



on the right hand-side of Eq. ©. Therefore, the QCD 



7 The spectral function in Eq. 1231 is supposed to capture 

the known perturbative contributions to d(q 2 ,m 2 r ). 



effective charge itself has to satisfy the integral represen- 
tation 



a(q 2 ,ml) 



g(a) 



da 



4m? 



(24) 



in this case as well. Otherwise, one would encounter 
a contradiction between the dispersion relation for the 
Adler D function of Eq. and its approximation given 
in Eq. ©. Besides, the condition l|24|) enforces the valid- 
ity of relations (0) and JHJ for the case of the nonvanishing 
pion mass. 

In general, there are several models for the invariant 
charge within the analytic approach to QCD (see Sec- 
tion [ffi] for the details). This is so by virtue of the fact 
that the behavior of the strong running coupling a s {q 2 ) 
at the ultraviolet asymptotic, which is known from per- 
turbation theory, together with the analyticity require- 
ment of the form of Eq. © or Eq. (|24|l . is not enough 
to uniquely determine the relevant spectral density p(a). 
The model for the analytic invariant charge 0, Il2l | has 
proved to be successful in the description of the strong 
interaction processes of both perturbative and intrinsi- 
cally nonperturbative nature |14l 1271 ] . We shall therefore 
adopt the spectral density of Eq. (|17|l in what follows. 

Thus, one arrives at the following integral representa- 
tion for the massive analytic invariant charge (see also 

Refs. mm) 



> m * ) = ToL 



a + z 



da, 



= X2> ( 2 5) 



where p [t '{a) denotes the ^-loop spectral density of 





Eq. (|17|l and x — Am 2 /A 2 . It is worth noting from 
the very beginning that the nonvanishing mass of the 
7r meson drastically affects the low energy behavior of 
this strong running coupling. Indeed, instead of the in- 
frared enhancement in the massless case of Eq. <|14|) . one 
has here the infrared finite limiting value for the massive 
invariant charge in Eq. I|25|l . 



£ f>(-)-, 
Po J y cr 



(26) 



which depends on the value of the pion mass. At the 
ultraviolet asymptotic, where the nonperturbative con- 
tributions are negligible, the result of Eq. t2*5)) tends to 
the perturbative running coupling ai (q 2 ). 



ailV,m^ a (, V) + 



rA 2 A 2 



4m, 



q 2 ln(q 2 /A 2 ) ' q 2 



(27) 



In this equation the limits q 2 A 2 and q 2 Am 2 , are 
assumed. In particular, the one-loop effective charge of 




<1 ReC 



FIG. 3: The integration contour in Eq. (0 for the case of 
the massive pion. The physical cut of the effective charge 
a(— (,mj) (see Eq. 1241 1 is shown along the positive semiaxis 
of real £. 



Eq. JUJ) reads for q 2 > 4m 2 



4 I 

(3 z 
4 I 
Po z 



4tt 1 

po z\nz 

7r / In x 
— + (x + 1) arctan 

2 V 7T 



e y arctan 



-) dy. (28) 



It is worthwhile to mention also that in the limit of mass- 
less pion TOjr = the effective charge (|25|1 coincides with 
the running coupling of Eq. (|I4|I . 

In order to handle the strong interaction processes in- 
volving the timelike kinematic variable one first has to re- 
late the experimental data with the perturbative results 
(see Section |nj. For practical purposes it is convenient 
to employ here the extension of the spacelike running 
coupling to the timelike domain given by Eq. J7J. The 
analytic properties in the q 2 variable of the QCD invari- 
ant charge a(q 2 ) are different for the massless © and 
massive (|2^)l cases (see Figures ^ and respectively). 
Thus, the continuation (J7J| of the massive strong run- 
ning coup ling 112511 to the timelike region results in (see 
Refs. Ha|4<| also) 



(s,m 2 ) 



PO Jw cr A 2 

(29) 

where s = —q 2 > 0, 9{x) stands for the Heaviside step- 
function (see, e.g., Rcf. 50]), p^\a) is the ^-loop spectral 
density defined in Eq. Q17[l. and x — Am 2 /A 2 . 

Let us address now the basic features of the running 
coupling in Eq. (|29p. First of all, it is very interesting 
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spacelike 
region 




FIG. 4: The one-loop massive analytic effective charge in the 
spacelike and timelike domains (Eqs. 1251 and 129L respec- 
tively). The values of parameters are: A = 623 MeV, ri{ — 2 
active quarks. 



to note here that the effective charges of Eq. 12ol) and 
Eq. ((23 have a common finite value in the infrared limit 
\q 2 \ — > 0, given by Eq. (SSjl. Second, the timelike mas- 
sive effective coupling of Eq. (|29[1 has the "plateau-like" 
behavior in the deep infrared domain: 



(30) 



Besides, for y/s > 2m w there is no difference between 
the massless and massive timelike running couplings of 
Eq. I|20l) and Eq. 1291) , respectively, since the mass of the 
7T meson affects the timelike effective charge ()29|l only in 
the region yfs < 2m^, where it does not run (see Fig- 
ure Therefore, in the ultraviolet asymptotic s — > oo, 
the expansion (|21|l . which accounts for the 7r 2 -terms, also 
holds for the running coupling of Eq. <|29|) . The hypoth- 
esis due to Schwinger |46|, H3] concerning the proportion- 
ality between the (3 function and the relevant spectral 
density holds for the massive timelike effective charge of 
Eq. I|29|) as well. Apparently, in the limit of vanishing 
pion mass — > the results of this section reproduce 
the massless case described in Section Hill 



It is worth noting here that some other models for 
the QCD effective charge also display a plateau similar 
to (|30f) in the infrared domain. In particular, the afore- 
mentioned optimized perturbation theory method [l6| 
predicts the stagnation of the timelike effective coupling 
in the region yfs < 300 MeV, in striking coincidence with 
the result obtained in Eq. i|5U|) . Moreover, the so-called 
"H-model" , with a similar freezing of the effective charge 
to a constant value in the infrared domain, has proved to 
be useful in studying of the dynamical chiral symmetry 
breaking (see, e.g., Ref. [HI). 



V. INCLUSIVE r LEPTON DECAY 

In order to draw a quantitative conclusion on the low 
energy behavior of a model for the strong running cou- 
pling, one needs the relevant estimation of the QCD scale 
parameter A. The latter can be extracted, for exam- 
ple, from the experimental data on the strong interaction 
processes. Among them, the measurement of the r de- 
cay width is most suitable for our purposes, since these 
data are fairly precise, and this process probes the in- 
frared hadron dynamics at energies below the r lepton 
mass < \/s < M T . Let us turn now to the study of 
this hadron process, restricting ourselves to the one-loop 
level at this stage. 

The experimentally measurable quantity here is the 
inclusive semilcptonic branching ratio 



R T 



r(r — > hadrons v T ) 
T(t- — > e~ v e v T ) 



(31) 



One can split this ratio into three parts, namely R T — 
R T . V + R t ,a + Rt.s- The terms R T . V and R t ,a account for 
the contributions to Eq. (|3*T|l of the decay modes with the 
light quarks only, and they correspond to the vector (V) 
and axial- vector (A) quark currents, respectively. The 
accuracy of the experimental measurement of these terms 
is several times higher than that of the strange width 
ratio R T ,s, which accounts for the contribution to Eq. I|31|) 
of the decay modes with the s quark. Thus, let us proceed 
with the nonstrange part of the ratio R T (|31|l associated 
with the vector quark currents 

Rt,v = ~y \V u d\ 2 Sew (l + S' EW + S QCD ) , (32) 

see Refs. [3^, [S^, E3 for detailed discussion of this is- 
sue. In Eq. (|32|l N c — 3 is the number of colors, 
|Kd| = 0.9738±0.0005 denotes the Cabibbo-Kobayashi- 
Maskawa matrix element S EW = 1.0194±0. 0050 and 
5' EVV — 0.0010 are the electroweak corrections [53,[55j, and 
(5 QCD stands for the strong correction. The recent mea- 
surements of the ratio 1)32(1 gave i? r ,v = 1.775 ± 0.017 
(ALEPH Collaboration, Ref. H3) and i? T , v = 1.764 ± 
0.016 (OPAL Collaboration, Ref. Assuming that 

these data have equal statistical weights, one arrives at 
the averaged value 



R TtV = 1.769 ±0.017. 



(33) 



In the framework of the approach in hand the strong 
correction in Eq. I|32|l at the one-loop level takes the form 



1 



M*J 



1 



2 W" 



ds 



M 2 ' 
(34) 

see, e.g, papers 0, [3t| and references therein. In 
Eq. I|34|) c^ 1 ) (s) is the one-loop strong running coupling 
in the timelike region, and M T = (1776.99to;2e) MeV 
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denotes the mass of the r lepton [5^|. As was shown in 
previous sections, the mass of the it meson entering the 
dispersion relation for the Adler D function (0 affects 
the low energy behavior of the QCD effective charge a(s). 
Consequently, handling the experimental data on the in- 
clusive r lepton decay is different for the cases of massless 
and massive pion. In order to demonstrate how the es- 
timation of the QCD scale parameter A is affected by 
the nonvanishing mass of the 7r meson, let us study both 
instances. 

For the limit of massless pion = 0, the one-loop 
strong correction S QCD to the R TiV ratio is given by 
Eq. J2U, with a^(s) being the one-loop massless effec- 
tive charge of Eq. I|20[l . Although the latter possesses 
the enhancement at s — > (see Eq. I|22flh the resulting 
singularity is integrable. Then, it is useful to represent 
the QCD correction in a more convenient form 

Po Jo 

+ isi 1 n )(Af r 2 ), (35) 

where ain (s) is the running coupling (|20[1 . p^(a) de- 
notes the one-loop spectral density (|19|l . and the nota- 
tions £ = s/AI 2 and fi — M 2 /A 2 are used. For the 
experimental data given in Eq. I|33[) one gets the value 
A = (508 ±61) MeV for the QCD scale parameter. This 
estimation corresponds to rif — 2 active quarks, and its 
uncertainty is due to the errors in the values of i? TjV , 
I Kid |, <Sewj and M T . The relevant behavior of the mass- 
less analytic invariant charge in the spacelike and time- 
like regions (Eqs. i|16[l and l|2U|) . respectively) is shown in 
Figure E 

Let us proceed now to the case of the nonvanishing 
pion mass. Here, the one-loop QCD correction to the 
i? TjV ratio reads as 

9 r M T / \ 2 

x(l + 2^)aS( 8 ,mi)^,(36) 

where ain (s, to 2 ) is the one-loop massive analytic charge 
of Eq. jjgfy and to„ = (134.9766 ± 0.0006) MeV stands 
for the 7T° meson mass [54| . In general, in the framework 
of the analytic approach there is no need to involve the 
contour integration in Eq. (|34|l , since the effective charge 
a?(s), appearing in the integrand, contains no unphysical 
singularities in the region s > 0. In other words, the inte- 
gration in Eq. I|34l) can be performed in a straightforward 
way. Thus, one can cast the strong correction l|36|) into 
a convenient form 

SiUM^ml) = ± /V-2£ 2 + 2)p«(^ 

+ isW(M T 2 ,TO 2 ), (37) 



where £o = 4m 2 /Af 2 and the other notations have been 
explained above. For the experimental data ((33(1 the es- 
timation 8 A = (623 ±81) MeV has been obtained for 
nf = 2 active quarks. The uncertainty here is because of 
the errors of i? TjV , |V u( j|, S EW , M T , and m^. The corre- 
sponding infrared limiting value of the massive effective 

charge (231 is a 1} = 1.475 ± 0.170. The low energy be- 
havior of the analytic running coupling in the spacelike 
and timelike domains (Eqs. (|25|l and 1)29(1 . respectively) 
is presented in Figure 

Thus, in the framework of the approach in hand it 
proves to be important to take into account the mass 
of the 7r meson in processing the low energy QCD data. 
Specifically, the relative difference between the obtained 
estimations of the scale parameter A for the massive and 
massless cases is about 20 %. This is so by virtue of the 
fact that the contribution to the strong correction l|34|) 
of the effects due to the pion mass 

A^oCM 2 , to 2 ) = ± i f - 2£ 2 + 2) p«(£ M ) d£ 
Po Jo 

(38) 

turns out to be significant. At the same time, since the 
scales involved in the integral (|38|l are very low, for some 
models for the analytic running coupling the difference 
between the limits of massive and massless 7r meson may 
not be so sizable. For example, in the case of the Shirkov- 
Solovtsov model p5|. where the relevant spectral den- 
sity reads as pil\a) = l/(ln 2 <7 + 7r 2 ), the relative dif- 
ference between the values of the QCD scale parameter, 
extracted from the experimental data on the inclusive 
r lepton decay (|33|l . is about 1 %, but the obtained esti- 
mations appear to be rather large. Namely, at the one- 
loop level with rtf = 2 active quarks one gets the values 
A = (9651^2) MeV for the case of the massless tt meson, 
and A = (976i 28 3) MeV for the nonvanishing pion mass. 



VI. APPLICABILITY TO THE CHIRAL 
SYMMETRY BREAKING 

Based on the study of gauge invariant Schwinger- 
Dyson equations, Cornwall proposed a long time ago that 
the self-interactions of gluons give rise to a dynamical 
gluon mass, while preserving at the same time the local 
gauge symmetry of the theory |2£j . This gluon "mass" is 
not a directly measurable quantity, but has to be related 
with other physical quantities, such as the glueball spec- 
trum, the energy needed to pop two gluons out of the 
vacuum, the QCD string tension, or the QCD vacuum 
energy (see paper jH^I and references therein). 



It is worthwhile to note here that the one-loop perturbative anal- 
ysis of the strong correction in Eq. 1321 (see, e.g., Ref. lo.'J i gives 
the value of the QCD scale parameter A = (690 ± 57) MeV for 
two active quarks. 
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FIG. 5: Graphical representation of the gap equation 14211 . 



One of the main phenomcnological implications of this 
analysis is that the presence of the gluon mass m g satu- 
rates 9 the running of the strong coupling at low energies. 
Namely, instead of increasing indefinitely in the infrared, 
as perturbation theory predicts, it "freezes" at a finite 
value, determined by the gluon mass. In particular, the 
nonperturbative effective coupling obtained in Ref. [2jJ 
is given by 



a c (q 2 ) 



4n 



1 



/3 m[z + 4M 2 ( g 2 )/A 2 



Si 

A 2 ' 



(39) 



where M g (q 2 ) denotes the dynamical gluon mass 



ln(z + 4m 2 /A 2 ) 
ln(4m 2 /A 2 ) 



(40) 



Here the nontrivial dependence of the dynamically gen- 
erated gluon mass (|40|l on the momentum q 2 is crucial 
for the renormalizability of the theory. The running cou- 
pling l|39|l has the infrared finite limiting value a c (0) = 
4-7T [/?o ln(4m 2 /A 2 )] . It is worth noting that the above 
equation makes sense only for the gluon mass satisfying 
m g > A/2. For a typical values of m g = 500 MeV and 
A = 300 MeV, one obtains for the case of the pure gluo- 
dynamics (rif = 0) an estimation a c . (0) — 0.5. An inde- 
pendent analysis presented in Ref. [60J yields a maximum 
allowed value for a c (0) of about 0.6. The incorporation 
of fermions into the effective charge [6l| does not change 
the picture qualitatively (at least for the quark masses of 
the order of A), resulting in an approximate expression 



a cp (q 2 ) 



4?r 



11 ln(z + Xg ) - 2n t \n{z + x,)/3 ' 



(41) 



In this equation Xg = 4to 2 /A 2 , Xq = 4m 2 ,/A 2 , a light 



quark constituent mass is m q = 350 MeV ,54], and 
m g = (500 ± 100) MeV stands for the gluon mass. The 
effective coupling of Eq. (|39|) was the focal point of ex- 
tensive scrutiny, and has been demonstrated to furnish 
an unified description of a wide variety of the low energy 
QCD data 

In general, an important unresolved question in this 
context is the incorporation of the QCD effective charge 



9 Another discussion of the impact of the gluon mass on the in- 
frared beha vior of the strong running coupling can be found, e.g., 
in Ref. IH3 . 




FIG. 6: A typical dependence of E(0) on the infrared limiting 
value of the QCD effective charge a(0), S~ 1 (p) = A(p)$ + 
E(p). 



into the standard Schwinger-Dyson equation governing 
the dynamics of the quark propagator S(p) 



S~ 1 (p)^S ( ; 1 (p)-g' 



d 4 k 



ltl SY v ^, (42) 



see Figure |S] also. In particular, since QCD is not a 
fixed point theory, the usual QED-inspired gap equation 
must be modified, in order to incorporate the running 
charge and asymptotic freedom. The usual way of ac- 
complishing this eventually reduces to the replacement 
1/fc 2 — * a(k 2 )/k 2 in the corresponding kernel of the gap 
equation, where a(k 2 ) is the QCD running coupling. The 
inclusion of a(k 2 ) is essential for arriving at an integral 
equation for S(p) which is well-behaved in the ultravi- 
olet. Indeed, the additional logarithm in the denomi- 
nator of the kernel due to the running coupling a(k 2 ) 
improves the convergence of the integral. However, since 
the perturbative form of a(k 2 ) diverges at low energies as 
l/ln(fc 2 /A 2 ) when k 2 — > A 2 , some form of the infrared 
regularization for the invariant charge a(k 2 ) is needed, 
whose details depend on the specific assumptions one is 
making regarding the nonperturbative hadron dynamics. 
At this point the issue of the critical coupling makes its 
appearance. Specifically, as is well-known, there is a crit- 
ical infrared limiting value of the running coupling, to be 
denoted by a cr , below which there are no nontrivial so- 
lutions to the resulting gap equation, i.e., there is no 
chiral symmetry breaking, see Figure Thus, the in- 
variant charge a(k 2 ) employed within the gap equation 
must be such that (i) it gives rise to a nonsingular an- 
swer, (ii) it reaches large enough values at k 2 — + in 
order to overcome a cr , and (iii) it does not contradict 
existing low-energy experimental results. 

The incorporation of the effective charge of Eq. I|39|) 
into a gap equation has been studied for the first time 
in Ref. |63j . There it was concluded that chiral symme- 
try breaking solutions for £(p) could be obtained only 
for unnaturally small values of the gluon mass, namely 
m g / A ~ 0.8. This is so because the typical value of a CI 
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0.0 



2 4 6 q, GeV 10 

FIG. 7: Comparison of the massive analytic running cou- 
pling (12ol (solid curves) with the effective charge l|41[l (dashed 
curves). The values of parameters are: rif = 2 active quarks, 
A = 704 MeV (a), A = 542 MeV (b), gluon mass m g = 
400 MeV and A = 350 MeV (c), gluon mass m g = 600 MeV 
and A = 150 MeV (d). 



found in the standard treatment of the gap equation 10 is 
cv cr ~ 1.2 (see, e.g., Ref. JtHj ), which is what the expres- 
sion for a c (0) yields for the above ratio of m s /A. This 
issue was further investigated in Ref. |6lj | , where a system 
of coupled gap and vertex equations was considered. The 
upshot of this study was that no consistent solutions to 
the system of integral equations could be found, due to 
the fact that the allowed values for a(0), dictated by the 
vertex equation, were significantly lower than ev cr , i.e., 
not large enough to trigger chiral symmetry breaking. A 
similar analysis was presented in Ref. |65| , together with 
several other models for the nonperturbative QCD run- 
ning coupling 66] . 

In what follows we will suggest a possible resolution 
of this problem, inspired by the infrared behavior of the 
massive analytic invariant charge of Eq. 1251) . The basic 
observation is captured in Figured Namely, the effec- 
tive charge with a gluon mass (dashed curves) and the 
analytic charge (solid curves) coincide for a large range 
of momenta, and they only begin to differ appreciably in 
the deep infrared domain k 2 < A 2 . In this region the 
analytic charge 125(1 rises abruptly, almost doubling its 
size between k 2 = A 2 and k 2 — 0, whereas the running 
coupling (|41|l in the same momentum interval remains 
essentially fixed to a value 11 of about 0.6. A possible 



10 The exact value of a cr depends on the number of active flavors 
as well as on the various approximations employed in deriving 
the gap equations, such as the choice of gauge, or the inclusion of 
gauge-technique inspired Ansatze for the quark-gluon vertex, but 
these issues do not alter significantly our qualitative discussion. 

11 The precise numerical values of m g /A and m n /A do not change 



picture that emerges from this observation is the follow- 
ing. It may be that the concept of the dynamically gener- 
ated gluon mass fails to capture all the relevant dynam- 
ics in the very deep infrared, where confinement or other 
nonperturbative effects make their appearance. At that 
point it could be preferable to switch to a description 
in terms of the analytic charge l|25|) . which (i) coincides 
with that of Cornwall in the region where the latter fur- 
nishes a successful description of data, and (ii) since it 
overcomes the critical value a cl ., offers the possibility of 
accounting for chiral symmetry breaking at the level of 
gap equations. 



VII. CONCLUSIONS 

In this paper the effects due to the mass of the 7r me- 
son are incorporated into the analytic approach to QCD. 
The nonvanishing pion mass gives rise to an infrared fi- 
nite limiting value for the QCD effective charge. Besides, 
the latter acquires the plateau-like behavior in the deep 
infrared domain of the timelike region < y/s < Ira^. It 
is of a particular interest to note that such stagnation is 
also predicted by a number of phenomenological models 
for the strong running coupling. The developed analytic 
effective charge is applied to processing the experimen- 
tal data on the inclusive r lepton decay. The effects due 
to the pion mass play a substantial role here, affecting 
the estimation of the QCD scale parameter A. A quan- 
titative conclusion on the applicability of the obtained 
massive running coupling to the study of chiral symme- 
try breaking is drawn. 

It would be interesting to further scrutinize the devel- 
oped approach. First of all, it is of particular relevance to 
include the higher order perturbative corrections in the 
study of the experimental data on the inclusive r lepton 
decay. Moreover, it might also be important to incor- 
porate the nonperturbative terms, arisin g fr om the op- 
erator product expansion (see also Ref. l67l ) and from 
the so-called nonlocal chiral quark model [63, into the 
Adler D function. In addition, a detailed study of the 
gap equation, with the analytic charged plugged into it, 
is needed in order to verify if indeed one encounters non- 
trivial solutions, whose size is phenomenologically rele- 
vant. Specifically, one should check by making use of, 
e.g., the Pagels-Stokar method 69], whether the solu- 
tions obtained for can reproduce a reasonable value 
of the pion-decay constant f n . It is also interesting to ap- 
ply the developed model to the study of the pion electro- 
magnetic form factor F v (q) (see paper [TtJ and references 
therein). At the same time, a crucial point to explore is 
whether the massive analytic effective charge satisfies a 
variety of phenomenological constraints, imposed by the 



qualitatively this picture, as long as the two scales are well sep- 
arated. 
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low energy experimental data on the infrared behavior of 
the QCD running coupling, see, e.g., papers [62,[71| and 
references therein. 
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